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1. Introduction 

The nature and origin of fermionic pairing was first elucidated in the pioneering 
work by Bardeen, Cooper, and Schrieffer, widely known as the BCS theory pQ. The 
attractive pairwise interaction between electrons with opposite spin, albeit extremely 
weak, can give rise to an instability in normal electron gas towards the formation of 
zero-momentum Cooper pairs near Fermi surface, and because of pair condensation, 
ordering of conduction electrons emerge naturally. When subjected to an external 
Zeeman field, the population balance between electrons with different spins is broken. 
As a consequence, not all electrons can find a partner to pair up with. If spin-population 
imbalance is large enough, the pairing of fermions has to occur at finite center-of- 
mass momentum with deformed Fermi surface state [2]. This exotic possibility of 
inhomogeneous superfluid was first predicted by Fulde and Ferrell (FF) [3j, and Larkin 
and Ovchinnikov (LO) [4] a little later. FF refers to an order parameter with plane- 
wave form A(r) = A e zqr , which spontaneously breaks time-reversal symmetry; while 
LO considers the superfluid with a standing- wave order parameter A(r) = A cos(q • r), 
which explicitly breaks translational symmetry. Both phases have puzzled the solid-state 
community for decades in terms of unambiguous experimental evidence to prove their 
existence. Moreover, the FFLO state is also of interest in quantum chromodynamics at 
low temperature and high density, where the property of asymptotic freedom may favor 
color superconductivity [5]. 

In recent years, due to their exquisite controllability, ultracold atoms have emerged 
as an ideal platform to simulate real materials. Adjustable interaction and high degrees 
of control over spin-populations have enabled one with the feasibility to explore this 
long sought phase. Tremendous theoretical and experimental efforts have been put into 
optimizing the best detectable parameter regime of FFLO phase. The most promising 
route is now believed to probe the one-dimensional (ID) spin-imbalanced Fermi gas 
[3 [3 EJ E] , where indirect evidence of FFLO phase has been found in a recent experiment 
[TO] . However, for 3D Fermi gas, the FFLO phase is not favored [TT1 [T2l [T3] . 

Over the past few years, another milestone achievement in cold atom research is 
the realization of artificial spin-orbit (SO) coupling, first in bosonic systems [14J and 
later in fermioinic ones [T5l [16] . By tailoring the laser fields that generate the SO 
coupling, various coupling schemes can be realized in principle. It has been realized very 
recently that, in a Fermi gas, the interplay between the SO coupling and an effective 
Zeeman field may lead to distortion of single-particle dispersion as well as the Fermi 
surface, in such a way that finite-momentum dimer state and/or Cooper pairs will 
be favored [T71 [HJ [THJ, [20j [21] . In this work, we provide a unified treatment of both 
two-body and the many-body physics for a Fermi gas subjected to an isotropic three- 
dimensional SO coupling (3DSOC) and an effective Zeeman field. The generation of 
such 3DSOC has been recently proposed by optically dressing four internal atomic states 
with a tetrahedral geometry [22J, [23]. This version of the SO coupling is less explored 
and unfamiliar to condensed matter community, where 2D Rashba and Dresshauls SO 
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couplings are studied extensively. One important advantage of 3DSOC over lower- 
dimensional SO interaction is that it provides the greatest enhancement of fermionic 
pairing. Furthermore, due to its isotropic nature, mathematical simplicity is ensured. 

The main findings of our work are: i) Under arbitrarily weak Zeeman field, zero- 
momentum dimer state and conventional BCS superfluid phase are no longer stable, 
ii) For many-body system, FF state is inherently robust, and ultimately connects to 
normal phase in a smooth manner as Zeeman field strength is increased, cf. Fig. [3} 
Moreover, this type of exotic superfluid has different origin in comparison with the 
previously studied FFLO state. In the absence of the SO coupling, individual particle 
number with different spins is conserved, hence the imbalance induced finite-momentum 
pairing has parity symmetry between q and — q, which should be called LO phase 
by definition; on the other hand, in the presence of the SO coupling, Zeeman field 
breaks time-reversal symmetry explicitly and cause the single-particle dispersion to be 
asymmetric, which underlines the idea of finite-momentum dimer bound state [17J and 
the FF pairing instability. The center-of-mass momentum of the Cooper pair can be 
as large as the Fermi momentum. This result should be very encouraging for future 
experimental exploration. 

The rest of the paper is organized as follows. In Sec. [2| we formulate the physical 
model and introduce the Functional Path Integral technique. We apply this general 
approach to the system with 3DSOC and discuss our results on two-body physics in 
Sec. [3] and on many-body physics in Sec. [4| And finally we conclude in Sec. [5| 

2. Physical model and general formalism 

In this section, we first present the model Hamiltonian under study, and then introduce 
the widely used functional path integral approach. Using this approach, we can 
discuss both the two-body physics and the many-body physics at both zero and finite 
temperatures in a unified way. 

2.1. Model Hamiltonian 

We start by formulating the Hamiltonian for a non-interacting homogeneous Fermi gas 



where £k = h 2 K 2 /{2m) — \i and ip = [^(r),^(r)] T is the fermionic annihilation 
field operator. We have defined SO coupling strength vector v = (y x ,v y ,v z ), and 
the Zeeman field vector A = (A x , A y , A z ). a = [o x ^a v ^a z ) are Pauli matrices acting 
on the atomic (pseudo-)spin degrees of freedom. This description is a general model 
valid for various SO coupling schemes. The single particle spectrum is given by 
E 1 (k) = + 7\/J2i( v iki + ^i) 2 with 7 = ±1 denoting the two helicity branches. 
The experimentally realized [TU [151 HH] equal weight Rashba-Dresselhaus SO coupling 



in 3D: 




(1) 
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takes the form of Eq. ([!]) with v = (0, 0, h 2 k r /m) and A = (£2/2,0,5) where k r is the 
laser recoil momentum, Q the Raman laser coupling strength and S the two-photon 
detuning. The Rashba SO coupling (25J [26] can be recognized with v = (v x , v y , 0) 
and v x = v y = h 2 k r /m. In our work, we will focus on the 3DSOC [22l [23] with 
v x = v y = v z = v. For this case, due to the isotropic nature of the SO coupling term, 
the direction of the Zeeman field is irrelevant and we shall choose it to be along the 
z-axis, and hence A = (0, 0, h). 

It is important to note that the Zeeman field does induce an asymmetry in the 
single-particle dispersion relation. To illustrate this, we consider a filled Fermi sea with 
simple topology (cf. [21]) at zero temperature. In Fig. [TJ we plot Fermi surface without 
and with Zeeman field. In the absence of the Zeeman field, the Fermi surfaces for both 
helicity branches are represented by spheres centered at zero momentum, as shown in 
Fig. [Ifa). When we turn on the Zeeman field, both Fermi surfaces are distorted and 
no longer possess reflection symmetry about the k z = plane, as can be clearly seen 
in Fig. [jjb). In this perspective, the ground state is associated with nonzero total 
momentum along the fc^-axis. 

Next we consider the attractive s-wave contact interaction between un-like spins 
which, in terms of the creation and annihilation field operators for the original spin 
states, is represented by 

Hint = U f dr^{(r)^(r)^(r)^ t (r) (2) 

where Uq is the bare coupling strength to be renormalized using the s-wave scattering 
length a s . In this work, we constrain our attention to the experimentally exploited 
broad Feshbach resonances, which is well captured by the single-channel Hamiltonian 
prescribed above. 

2.2. Functional Path Integral Formalism 

In this section, we describe the functional path integral technique [271 EHl EH] and 
start from the partition function Z — J V[ip (r, r) , $ (r, r)] exp {— £ [ip (r, r) , j/j (r, r)] } 
where the action 

S ty, $\ = f P dr fdrV ^ (r, r) d T ^ (r, r) + U (^, (3) 
Jo [J a 

is written as an integral over imaginary time r. Here /3 = 1/(&bT) is the inverse 
temperature and % , / 0) is obtained by replacing field operators ^ and ip with 
grassmann variables ijj and ^, respectively. We can integrate out the quartic interaction 
term using the Hubbard-Stratonovich transformation [29j, from which the pairing field 
A (r, r) is defined. If we assume the mean-field order parameter to be of FF-type 
A = A e 2q r , and further integrate out the fermionic fields, we arrive at an effective 
action as 

5U = f drf dr(-M) - ^log[-g A '] + PY. — + 2~" +< "' - < 4 > 
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Figure 1. (a) Fermi surfaces (a cut in the k y = plane) in the absence of the Zeeman 
field. The two concentric Fermi surfaces are spherically symmetric. The inner blue 
sphere represents the Fermi surface of the + helicity branch, while the outer yellow 
sphere of the — helicity branch, (b) Fermi surfaces in the presence of the Zeeman field 
along the z-axis: both Fermi surfaces are deformed in such a way that the cylindrical 
symmetry about the /c^-axis is still preserved, but the relfection symmetry about the 
k z = plane is broken. 



- Ck+q/2 /-+ 

-iA a y 



lOJ ri 



iA a, 

■6 



k-q/2 



f- 



(5) 



where f± = (vi(ki ± ^) ± A$) cr^. In the second term of Eq. ([4]), the trace is to be 
taken over the Nambu spinor space $ (r,r) = , V^] T , the real coordinate space 

and imaginary time. The last term in Eq. ^ comes from interchanging fermionic fields 
^ and ^with ^ and ^ and the corresponding equal-time limiting procedure [29J. 
From Eq. Q, we can further sum over Matsubara frequencies to arrive at the grand 
thermodynamic potential 



CI 1 \A\ 2 m 1 ( gk+q/2 + g-k+q/2 1 A , , 

^ s k \ a=l 

+^-^E"»0+«PH«>)) (6) 

where we have regularized the bare interaction strength C/ in terms of the s-wave 

scattering length a s by ^ = j^p^ ^Sk2^- are the quasi-particle energy 

dispersion, which are just the four eigenvalues obtained by solving det[^ 1 (k, Eg)] = 0. 
In our case, Eg are too complicated to be analytically presented here. 

In the following sections, we restrain our attention to the isotropic 3DSOC with 
Zeeman field along the z-axis, i.e., A = (0, 0, h) ) and use the ansatz for FF order 
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Figure 2. Finite momentum dimer bound state solution for 3DSOC. The coloring 
shows the magnitude of #2b, varying with the SO coupling strength v and the Zeeman 
field strength h\ the inset shows bound state energy E^ as a function of q/kp (along 
the z-axis) for different h, from top to bottom h/Ep = 0,0.1,0.2. We fix scattering 
length as l/(fc j pa s ) = — 1 



parameter A(r) = Ae iqz . 

3. Results on two-body problem 

Following the path integral approach, we can characterize two-body properties at low- 
energy sector by inverse vertex function, which we refer to [30j for more details. We 
found consistent results reported in our previous paper [17] which are obtained by solving 
the two-body Schrodinger equation. For a bound-state with total momentum q, the 
corresponding energy £" q is obtained by solving the following equation: 



rn 



Wa s = V ^ 1 c _ ( gL k 2 , k 2\ + 2e k f (7) 

k [ ^ ^k, q \e* ftl -(2h+vq)* K± + Kz) J 

where = — eq +k — ea_ k and e k = h 2 k 2 /(2m). For a given set of parameters /i, 
v and a s , we can numerically obtain the eigenenergy of the dimer bound state £" q as a 
function of q. The momentum q at which i? q reaches the minimum labels the dimer 
state with lowest energy. The binding energy is defined as = 2£ , min — i? qo , where E m [ n 
is the ground state energy of single particle spectrum ^"(k). Only when > can we 
consider the dimer as a true two-body bound state. Otherwise, its energy lies in the 
single particle continuum. For the convenience of further comparison with many-body 
state, here we take the laser recoil momentum k r: which determines the SO coupling 
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strength, to be equal to Fermi momentum fcp, which is determined by typical atomic 
density in experiments [HJ HHJ [IS] . 

For a Zeeman field along the z-axis, we have q = q 2 bZ- Following above- 
mentioned protocol, we plot q 2 b as a function of coupling strength and h in Fig. [2| 
As one would expect, Zeeman field tends to destroy two-body bound state; whereas SO 
coupling enhances its formation. The competition between these two outlines the critical 
boundary value, beyond which 65 becomes negative and no stable bound state can be 
found. With increasing h : the minimum of F q deviates away from zero momentum to 
some finite value. As long as Zeeman field is non-zero, the lowest-energy bound state 
would occur at finite center-of-mass momentum q 2 b- Our calculation shows that the 
magnitude of q 2 b can be as high as 0.2fcp. 



4. Results on many-body problem 

Motivated by the two-body results, one naturally attempts to explore the direct analog 
for the many-body system, which we study in this section. 

We take a canonical ensemble approach by considering a homogeneous system with 
fixed particle number N and volume V, and hence the density n = N/V = fc^/(37r 2 ). 
The important quantity that determines the mean field phase diagram shall be the free 
energy, also known as Landau potential, defined as F = Q + fiN. At zero temperature, 
it coincides with the ground state energy. For a given set of parameters (including 
SO coupling strength v, Zeeman field strength h, interaction parameter l/(fcpa s ), and 
temperature T ), order parameter A, chemical potential /i, and the FF momentum 
q = q^F z should be determined self-consistently by stationary conditions 
OF OF OF 

We shall explicitly consider three types of phase: normal gas (A = 0, q — 0), BCS state 
(A ^ 0, q = 0), and FF state (A ^ 0, q FF ^ 0). 



4.1. Zero-temperature phase diagram on the BCS side 

We shall first focus on a relatively weak-interacting regime on the BCS side of crossover 
and take l/(fcpa s ) = —1. In this regime, we can easily justify the mean-field treatment 
at both zero and finite temperature, and furthermore the SO coupling effect wound be 
more pronounced j34j |35l [36J [37] . 

To get some insights first, in Fig. [3j we plot free energy as a function of h for a 
given SO coupling strength v = E F /kp. (We choose this relatively large SO strength, 
to avoid possible complications, e.g. Sarma phase [38], phase separation [39] etc.) It 
is very remarkable to notice that FF state is energetically favored for arbitrarily small 
h. For instance, at h — 0.02£' j p, the gain of energy over the BCS pairing phase is 
AF = F B cs — ^ff ~ 4.54588 x 10 -5 Ai?F. However, this energy gain quickly increases 
as h is increased. For example, at h = 0.28Ep : we have AF = 1.13728 x 10 -2 AFp 



Fulde-Ferrell Pairing instability in Spin-Orbit Coupled Fermi Gas 



8 




Figure 3. Zero temperature free energy comparison for fixed SO coupling strength 
v = Ep/kp and interaction strength l/kpa s = —1. FF superfluid phase is favored 
throughout the plotted parameter space. In the inset, we plot the BCS order parameter 
Abcs, the FF order parameter Aff (both in units of Ep), and the FF momentum q 
(in units of kp) as functions of h. 



which is more than two orders of magnitude larger and represents a very large energy 
value on the BCS side of Feshbach resonance. Once again, the idea of favoring FF phase 
is backed by the picture of Fermi surface deformation (cf. Fig. [T]) and two-body bound 
state solutions [Eq. and Fig. [2]. When we further increase h, BCS superfluid is taken 
over by normal phase as the BCS order parameter drops to zero rather sharply (see the 
inset of Fig. [3]); on the other hand, FF state connects to normal phase very smoothly 
at a much larger value of h. 

The FF state here has different origin with the conventional FFLO states in the 
absence of the SO coupling j6j U\ E], in which case, for a given interaction strength 
and with increasing population imbalance (i.e., Zeeman field), one would expect that 
competitions among various quantum phase (BCS, Sarma, FFLO, and normal phases) 
could lead to both first- and second-order phase transitions. By contrast in the presence 
of SO coupling, especially 3DSOC, FF state dominates almost the entire phase diagram, 
as we map out in the v-h plane Fig. [|(a). The BCS phase only exists on the axis (i.e., in 
the absence of either the Zeeman field or the SO coupling). Normal phase and FF phase 
are connected by a smooth boundary, which we identify by setting a threshold value 
of energy difference |F F f — ^normal ~ 10~ 5 NE F . Note that close to the boundary, A FF 
also becomes exceedingly small. For illustration purposes, we schematically added two 
small regions near v = 0, the LO (green) and the phase separation (blue) regions, in the 
phase diagram. The boundaries of these two phases in the absence of the SO coupling 
(i.e., at v = 0), which are well studied, are obtained from previous results [SHI EHJ [40] . 
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Figure 4. (a) Zero-temperature phase diagram at l/kpa s = — 1 in the parameter 
space spanned by h and v. The FF phase is divided into gaped and gapless region by 
the green dashed line. BCS state only exists strictly on the axis marked by two red 
straight lines. The two blue lines indicate the smooth boundary between FF state and 
normal phase. Within the FF phase, the color scale indicates the momentum q FF . LO 
and phase separation (PS) regions are added schematically for illustration purpose, 
(b) Single-particle excitation gap AE, FF order parameter A FF and momentum q FF 
as functions of h. The SO coupling strengths are v = 0.2Ep/kp (red curves) and 
0.5Ep/kp (black curves), (c) Aff and q^F as functions of v for h = 0.28Ep. In all 
plots, the energy is in units of Ep, and momentum in units of kp. 



It has been shown that, with increasing SO coupling strength, both these phases are 
suppressed rather rapidly [40] . 

Furthermore, in Fig. [3](a), FF phase is divided into the gaped and the gapless 
regions by examining the single-particle excitation gap AE = mm{\E£\}, where Eg are 
quasi-particle dispersions introduced in Eq. ([6]). As shown in Fig. |i](b), AE decreases 
monotonically as a function of h and drops to zero at some critical value of h c which 
depends on the SO coupling strength v. The critical value h c is represented by the 
green dashed line in Fig. [l](a). At h C) both q FF and A F f exhibit kinks for relatively 
small SO coupling strength. These kinks get washed out quickly with increasing v (see, 
for instance, the inset of Fig. [3]). On the other hand, in the limit of v = 0, these kinks 
become true jumps signaling the first-order phase transition between the BCS phase 
and the FF phase region. In Fig. [4](c), we plot q FF and A FF as functions of v for a fixed 
h. We note that even though A F f increases monotonically as v : the FF momentum q FF 
shows non-monotonic behavior: it first increases and then decreases as v is increased 
from zero. 

It is instructive to make comparisons between the two-body results and the many- 
body results. To this end, we consider a cloud of degenerate Fermi gas typically realized 
in experiment, with density n = 10 12 cm -3 which defines kp and Ep. We compare the 
two-body dimer momentum g 2 b with the many-body FF pairing momentum q FF in Fig. [5] 
at two different values of SO coupling strengths. Note that the range of h values for 
which the two-body bound state exists is much smaller than that for the existence of 
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Figure 5. Momentum comparison of the FF state and two-body bound state at 




Figure 6. Nodal Fermi surface plots in momentum space. The closed surfaces are 
formed by momentum values at which the excitation gap AE vanishes. Both figures 
has the same interaction parameter l/k F a s — —1 as in phase diagram Fig. [3] and 
coupling strength v = O.^E^/hp, while h = 0.2E F for (a) and h = 0.35E F for (b). 

the FF state. For example, at v = 1.3E F /k F) two-body bound states only exist for 
h < 0.1E F ; while the FF state extends all the way up to about h w E F . As such, the 
largest q^F that can be achieved is much larger than the largest g 2 b- I n the region where 
both two-body bound state and the FF state exist, q^F and q 2 b are comparable, with 
the latter somewhat larger. The difference between them becomes smaller as the SO 
coupling strength increases, indicating that at large SO coupling strength, the many- 
body properties of the system are also dominated by the two-body physics. 

Before ending this subsection, we want to remark on the gapless FF state. For 
Zeeman field strength above the critical value h C) one or more quasi-particle energy 
E£ will vanish at certain values of momentum k. Such momenta form closed surfaces 
(nodal Fermi surface) in momentum space with cylindrical symmetry around the k z - 
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Figure 7. Zero temperature phase diagram with interaction strength l/kpa s = —2 
and l/hpGLs = for (a) and (b), respectively. The color scale represents <7ff/&f- (c) 
FF superfluid momentum q^F as a function of the interaction strength. For curves 
from top to bottom, h/Ep = 0.25, 0.2, 0.15, 0.1, 0.05, and the SO coupling strength 
is fixed to be v = Ep/hp. 

axis and reflection symmetry about the k z = plane. Hence such nodal Fermi surfaces 
always appear in pairs and may be measured using the technique of momentum-resolved 
radio-frequency spectroscopy. Two examples are illustrated in Fig. |6| 

4.2. Effects of interaction 

So far we have focused on the zero-temperature phase diagram of a weakly-interacting 
system. Now we briefly discuss the effects of interaction and finite temperature. In 
Fig. [7|(a) and (b), we present two zero-temperature phase diagrams in the h-v plane for 
l/kpa s = —2 and 0, respectively. They are qualitatively similar to the one presented in 
Fig.gfa) for l/k F a s = -1. As we move from the BCS limit towards unitarity, the region 
for normal phase shrinks and the FF superfluid remains dominant. Furthermore, the 
region for gaped FF phase increases quickly. At unitarity, the whole parameter space 
presented in Fig. [7|(b) are occupied by the gaped FF phase. On the other hand, for fixed 
h and the FF momentum q FF quickly decreases as we go from the BCS side to the 
BEC side, as shown in Fig. [7|(c). This result is consistent with the one obtained from 
the two-body study [T71 HE] . 

4.3. Effects of temperature 

Finally, we consider the effects of finite temperature. In Fig. |8](a), we plot the phase 
diagram in the parameter space spanned by h and T by taking l/(kpa s ) = —1 and 
v = Ep/kp. The FF superfluid phase dominates at small h and low T. There is a 
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Figure 8. (a) Finite temperature phase diagram at l/kpa s = — 1 and v = Ep/kp. 
The color scale indicates FF order parameter Aff- (b) The free energy difference 
between the two superfluid phases (BCS and FF) and the normal phase at temperature 
T = O.lTp. FF phase always has lowest free energy. 

second order transition towards normal phase as h and/or T increases. The BCS phase 
again only lives on the h = axis. In Fig. [8](b), we compare the free energies for all 
three phases at T = O.lTp and clearly show that the FF phase possesses the lowest free 
energy at any finite values of h as long as h is below a threshold at which the system 
turns normal. 

5. Conclusion 

In summary, we have studied spin-orbit coupled Fermi gas subjected to an effective 
Zeeman field. Based on the picture of Fermi surface deformation, the two-body 
calculations, and the mean-field many-body results, we conclude that the BCS state 
with zero-momentum Cooper pairs is not stable against Fulde-Ferrell superfluid pairing 
at any finite Zeeman field strength. The FF phase is robust against interaction and 
finite temperature, and the corresponding center-of-mass momentum of the Cooper pair 
can be comparable to the Fermi momentum. 
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